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In the effective mass approximation, electronic property in graphene can be characterized by the 
relativistic Dirac equation. Within such a continuum model we investigate the electronic transport 
through graphene waveguides formed by connecting multiple segments of armchair-edged graphene 
nanoribbons of different widths. By using appropriate wavefunction connection conditions at the 
junction interfaces, we generalize the conventional transfer matrix approach to formulate the linear 
conductance of the graphene waveguide in terms of the structure parameters and the incident elec- 
tron energy. In comparison with the tight-binding calculation, we find that the generalized transfer 
matrix method works well in calculating the conductance spectrum of a graphene waveguide even 
with a complicated structure and relatively large size. The calculated conductance spectrum indi- 
cates that the graphene waveguide exhibits a well-defined insulating band around the Dirac point, 
even though all the constituent ribbon segments are gapless. We attribute the occurrence of the in- 
sulating band to the antiresonance effect which is intimately associated with the edge states localized 
at the shoulder regions of the junctions. Furthermore, such an insulating band can be sensitively 
shifted by a gate voltage, which suggests a device application of the graphene waveguide as an 
electric nanoswitch. 

PACS numbers: 84.40.Az, 81.05.Uw, 73.23.-b, 72.10.-d 



I. INTRODUCTION 

The fabrication of graphene, a graphitic sheet of one-atom thickness, is the first experimental realization 
of truly two-dimensional crystal^ Such a carbon material presents many unusual electronic and transport 
properties, such as the half-integer quantum Hall effect the nonzero conductivity minimum at van- 
ishing carrier concent rat ion ^iL&£ i 1 o,n 1 i2 1 i3 subtle weak localization ) 14 : 15 ' 16 : 17 ' 18 : 19 and the refiectionless 
transmission of the carrier through an arbitrarily high barrierj 20 ' 21 ' 22 From the application point of view, 
graphene has very high mobility and shows am-bipolar behavior^ More importantly, the planar geometry 
of graphene makes it relatively easier to fabricate various integrated nano-circuits. Therefore, graphene is 
regarded as a perspective base for the post-silicon electronics. 

Motivated by possible device applications, electronic transport through various graphene nano-structures 
were extensively studied both experimentally and theoreticall y 24 ' 25 ' 26 i 27 ' 28 i 29 ' 30 i 31 ' 32 ' 33 ' 34 Among these struc- 
tures, graphene nanoribbon(GNR) is the basic element to carry the current flow. Band structure calculation 
indicates that the zigzag-edged GNR is always metallic while the armchair-edged ones are either metallic 
or semiconducting, depending on their widths Recent experimental work has confirmed the possibility of 
tunnelling the transport gaps of GNRs by changing their widths^ 6 - 

A graphene junction can be formed by interconnecting two semi-infinite GNRs with different widths. 
In such a graphene nanostructure, a traveling carrier is scattered by the junction interface. Thus, a fi- 
nite junction conductance appearsj 32 ' 33 Furthermore, when connecting multiple segments of distinct GNRs 
in cascade manner, one can build up graphene multiple junction(GMJ) structure. It was experimentally 
demonstrated that the ballistic transport in graphene can be retained over sub-micron scaled Therefore, 
the scattering of geometrical shape of a GMJ plays a dominant role in determining its electronic transport 
property. Recently, the conductance spectrum of a graphene single junction is studied in details* 3 ^ 3 - It was 
found that the presence of the lattice vacancy can efficiently enhance the junction conductance, 34 because 
that a vacancy makes the coupling between the electron states of the two GNRs at the junction interface 
stronger. In contrast, the investigation on electronic transport through GMJs and other graphene waveg- 
uides is still in its infancy. However, one can reasonably expect that the graphene waveguides possess many 
interesting electronic transport properties controlled by their geometrical shapes. For instance, in contrast 
to a graphene single junction, the GMJ has more junction interfaces. Thus more electron partial waves, 
including the incident wave and the multiple reflected waves from different interfaces, will take part in the 
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quantum interference. As a result, some transport features observed in graphene single junction can be en- 
hanced or weakened, depending on the details of the geometrical shape of the GMJ. Prior to the experimental 
realization to fabricate graphene waveguides with high edge-order, theoretical prediction on the geometrical 
shape dominated electronic transport properties of typical graphene waveguides is highly desirable. 

So far two kinds of theoretical approaches, the tight-binding method and the continuum model, have been 
frequently employed to study the electronic and transport properties of graphene bulk material and graphene 
nanostructure . 35 ' 36 ' 37 ' 38 The tight-binding method is suitable to treat the electronic state of graphene 
nanostructures with arbitrary shapes. Combined with the Landauer-Biittiker formula in discrete lattice 
representation ; 39 ' 40 this method is convenient to study the electronic transport through graphene nanos- 
tructures, with various scattering mechanisms included naturally. However, as the size of the graphene 
nanostructure under consideration increases, the tight-binding method runs into an embarrassment since 
it needs to treat some large matrices. Hence the calculation becomes rather time-consuming. In contrast, 
such a problem does not occur in the continuum model which is established based on the effective mass 
approximation! 41 ' 42 The continuum model has succeeded in describing the electronic and transport proper- 
ties of the bulk graphene. When applying it to the nanostructures with relatively large sizes and regular 
shapes, some appropriate boundary conditions are needed to solve the Dirac equation in this model. 

In the present work, we will investigate theoretically the electronic transport through the GMJ structures. 
Such a graphene waveguide usually has relatively large size. Thus, we will adopt the continuum model in our 
theoretical treatment. To study the transport properties we known that the transfer matrix approach is well 
developed to describe the electronic transport through various quantum waveguides made of conventional 
semiconductor materials! 43 ' 44 ' 45 ' 46 However, when we attempt to apply this approach to graphene waveguide, 
some generalization is necessary since the electron in graphene obeys the relativistic Dirac equation, rather 
than the schrodinger equation for the conventional semiconductor materials. This is just the subject of our 
present work. By using our generalized transfer matrix method, we calculate the conductance spectrum of 
the GMJs. We find that the GMJs always present a well-defined insulating band around the Dirac point, 
no matter whether the GNRs in the GMJ are metallic or semiconducting. We analyze that the origin of 
the appearance of the insulating band is the antiresonance effect produced by the edge states localized at 
the lateral zigzag edges of the junction interface. Furthermore, we also find the position of the insulating 
band can be sensitively adjusted by exerting a gate voltage under the GMJ. According to such a feature, we 
suggest that the GMJ structures can be considered as a device prototype of a nanoswitch. 

The rest of this paper is organized as follows: In SecJITl a self-contained theoretical framework to formulate 
the linear conductance of the GMJ structures is elucidated, though it is, in fact, a generalization of the con- 
ventional transfer matrix method. At first, we establish the relationship between the linear conductance and 
the scattering matrix in the GMJ. Then we develop the transfer matrix approach to work out the scattering 
matrix, which is realized by means of the wavefunction connection condition and boundary condition at 
the junction interfaces. In Sec lIIIl the numerical result on the linear conductance spectrum of some typical 
GMJs are shown. The antiresonance driven insulating band shown in the calculated conductance spectrum 
is discussed. Finally, the main results are summarized in Sec lIVl 



II. THEORY 



A. Scattering matrix and linear conductance 



The geometry and the honeycomb lattice structure of the GMJ under our consideration is illustrated in 
Fig.l. The relevant structure parameters are also explained in the same figure. Obviously, the basic elements 
consisting of such a GMJ are armchair-edged GNRs of different widths. The two semi-infinite ribbons at the 
left and right sides serve as two leads for electronic transmission. The other finite segments of GNRs are in 
between, constituting the device region where the electronic tunneling is scattered. It was previously shown 
that the continuum model can describe the electronic property of an infinitely long armchair-edged GNR 
very welL 30 ' 47 The electronic eigen-state in such a ribbon obeys the following Dirac equation subjected to a 
certain boundary condition. 
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where k± = — id x ± d y , 7 = ^/3tao/2 is the so-called Fermi velocity and t being the electron hopping energy 
between the neighbor lattice points. We choose the Dirac point as the energy reference point and in what 
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follows we use the units K = t = do = 1 ■ The electronic eigen-state in the above equation can be analytically 
solved, which is given by 
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And the corresponding eigen-energy is E = sj^/q? + fc?. 



±1 denotes the conduction and valence bands, 



respectively. In Eq.Q W denotes the width of the ribbon and its length in y direction has been normalized 
as unity. 9j is the angle between the wavevector k — (qj,kj) and the x axis. The transverse wavevector 
component qj takes some discrete values due to the quantum confinement. 
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When the ribbon width takes some particular values: W = 3m/2 with the integer m > 2, such a wavevector 
can be alternatively expressed as qj — jir/W. Thus, qj = is allowed, which implies a zero bandgap between 
the conduction and valence bands. In such a particular case the GNR presents a metallic behavior. 

In each segment of the GMJ as shown in Fig.l, the wavefunction of the electronic state with energy E can 
be expanded in terms of the eigen-states of the armchair-edged GNR with the same width. To be concrete, 
such a wavefunction expansion in Ath segment can be expressed as 
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where the wavevector kj(p m ) > and kj(p m ) = ~kj(p m ). Evidently, the first two terms in the right-hand 
side of the above equation stand for the propagating modes along the positive and negative directions, 
respectively. Apart from these propagating modes, we need to add the evanescent modes in the above 
expansion. The corresponding eigen wavefunction is given by 
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Of course, these evanescent modes are unphysical solutions in any uniform ribbon with infinite length. 
However, they are indispensable in the present GMJ structure to make the wavefunction continuity at the 
junction interfaces. Noting that in the left and right leads we should let (3 1 = a N+1 = to prevent the 
wavefunction from diverging. In Eq.(|U) the relation between the wavevectors and the incident electron energy 
E is 
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where V represents a constant potential to mimic a possible gate voltage exerted in this segment. For a 
given energy E the number of the propagating mode appearing in the expansion, denoted by J x , is finite. 
And it can be readily determined by the above dispersion relation. On the other hand, one can also see from 
the above dispersion relation that the number of evanescent modes is in principle infinite. We can actually 
add an appropriate number of the evanescent modes in the expansion, J x , to guarantee the convergence of 
the calculated results. Thus the total mode number in Ath unit is specified, J x = J x + J x . Such a quantity 
varies among the segments, depending on their respective widths. 

The central task in our theory is to find the scattering matrix which links the injected and the reflected 
wave amplitudes from the two leads to the device region. Namely 
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Where a 1 and b 1 (a JV and b N ) are column matrices consisting of the expansion coefficients in Eq.((4]), namely, 
the wave probability amplitudes of all the propagating modes in the left(right) lead. To determine the 
scattering matrix [S] one should make use of the wavefunction continuity at all junction interfaces. And 
the formal theory is referred to as the transfer matrix method. But before going into the details about 
the transfer matrix, we had better to present a brief statement about the relationship between the linear 
conductance and the scattering matrix, particularly suitable to the GMJ structure. 
First of all, we define a pseudo time reversal operator 



T = {a z ® I)C = 
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with C being the complex-conjugate operators. It is commutable with the Dirac Hamiltonian given in 
Eq.JI]). Such a commutation means that in any segment of the GMJ, besides ip , Tip x are also the possible 
wavefunction with the same eigen-energy. To be specific, Tip 1 and Tip N are allowable wavefunctions in the 
left and right leads. Owing to such an argument the following relation holds true, 



a 1 
b N 



[S] 



b 1 * 

N* 



a 



By comparing it with Eq.© we find that the scattering matrix satisfies 

[S\* = [S]-\ 

To calculate the probability currents in both leads, we obtain 
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and 
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In the above two equations the velocity operator is defined as v y — i [H, y] . The probability current conser- 
vation requires Jl = Jr, namely 
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If we define a new scattering matrix [S'] = [v] [S] [v] 1 in which the diagonal matrix [v] is defined as [v] 
diag[{sin2 {sin 5 6™}], we can derive from Eq. (fl4|) that the matrix [S'] is unitary, namely 



[s'][s'y = i 



(15) 



From the above relation and Eq. (|ll|) we can deduce an important conclusion about the matrix [S']. That is 
the self-transposability of it. 
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Now suppose that an incident electron comes from the Ith mode in the left lead. Then it has certain 
probability to penetrate into the jth mode of the right lead. The corresponding transmission coefficient is 
defined as the ratio of the incident and the penetrated probability current. By a straightforward calculation 
we find that such a quantity is associated with the scattering matrix [S'] . It is given by 
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Conversely, for an electron incident from the jth mode in the right lead to enter into the Ith mode in the 
left lead, we can work out the transmission coefficient in the same way. 
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The self-transposability of [S'] guarantees the symmetry of the transmission coefficients, i.e. 

T l ^ j =T j ^ l =T lj . (19) 

By virtue of such a relation, when a small bias voltage is established between the left and right leads, the 
net current can be compactly expressed as 

1=1 f J2 T ^U(E) ~ h{E)]dE, (20) 
3 jl 

where f L (E) and f R (E) are the Fermi distribution functions in both leads. Consequently, from the above 
expression we can readily extract an expression of the linear conductance. It is given by 

G=^Y?» = tY)W- (21) 

Jl jl 



B. scattering matrix and transfer matrix 

Next we need to work out the scattering matrix with the help of the connection condition of wavefunctions 
at the junction interfaces. To illustrate our derivation, we consider the wavefunction continuity at an 
arbitrary interface, say the Ath one, as shown in Fig. 1(b). It is expressed as 

^ x (x,y x )=ip x+1 (x,y x ), Q<x<W x+ i, (22) 

and 

^ x {x,y' x )=^ x+x {x,y' x ), 0<x<W x+1 , (23) 

where V'a(V'b) ^ s a two-component spinor, representing the wavefunction component belonging to A(B) 
atoms in Ath segment. Noting that y' x = y x — (\/3) _1 in the above equations, it arises from the spatial 
difference between the adjacent A and B atoms. However, such a trivial difference can be safely ignored since 
in the continuum model which is particularly appropriate to describe the large size structure, the envelop 
wavefunction in the preceding equations is slow-varying in lattice constant scale. At the shoulder region 
of the Ath junction interface, the wavefunction component corresponding to the truncated A atoms should 
vanish. It yields a boundary condition as below 

i> x (x, y x ) = 0, W x+1 <x<W x . (24) 

At other junction interfaces there are analogous wavefunction connection and boundary conditions. From 
these conditions we can derive a relation about the expansion coefficients in the adjacent segments. First of 
all, for convenience we adopt the following Dirac notations to denote the respective components of the two 
sublattices(A and B) in the GNR eigen-state. For A atom it is denoted as 

1 / pill x \ a 

\^s x lk x ) = ^== « e«.» (25) 



and the one for B atom is 



y/2Wx \ -e'^i x 



\XBs x lk x ) = -s x e ie "\XAs x lk x ). (26) 



Such a separated expression is equivalent to the eigen wavefunction given by Eq.([2]) except for the trivial 
change of the normalization constant. For the evanescent modes allowed in the junction structure, the eigen 
wavefunction given by Eq.((S]) can also be denoted in a separated form. 

\\As x mp x m ) = -jL= t e * 9m l | (27) 

and 

\XBs x mp x n ) = - Xm \XAs x mpi). (28) 
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Noting that to normalize these evanescent eigenmodes is not needed because it does not affect the final form 
of the scattering matrix. With these notations the wavefunction connection condition given by Eq. (|22p can 
be rewritten in an expanded form. 

J2[a X \XA s x lk x ) + b x \XA s x lk x )] +^[<|AA s A m^) + (3 x m \XAs x mp x m )] = 

I m 

$>, A+1 |A + lAs x+1 3 k x+1 ) + b x+1 \X + lAs x+1 jkf +1 )} + (29) 

3 

5> A+1 |A + lAs^np^) + /3 A+1 |A + lAs^np^ 1 )]. 

71 

In order to find a relation between these expansion coefficients, we act on the eigenkets in both sides of the 
above equation by the eigenbra (XAs x lk x \. Special attention should be paid to the issue that the integration 
over x involved in the inner product calculation should be restricted in the narrow region of the junction, 
namely, [0, WA+i] as labeled in Fig.l. However, considering the boundary condition given by Eq. (|24p , such 
integrations in the left-hand side of Eq. ([29f can be safely extended to total interface, i.e. [0, W\]. Following 
such a rule, we obtain 

a x + e -nk x y x b X = ^2[(\As x lk?\\ + lAs x+1 jk x+1 )a x+1 + (XAs x lk?\X + \As x+1 fk x+1 )b x+1 ] (30) 

3 

+ ]T[(AAsX A |A + \As x+1 np x+1 ) a x+1 + (XAs x lk x \X + \A S x+1 np x + 1 )l3 x+1 ]. 

n 

Instead of (XAs lk\ ; |, using the eigenbra (XAs x mp^ n \ to take the inner product we have 

e^al + /t = '^^[(XAs x mp x n \\ + \As x+1 ] k x+1 )a x+1 + (XA S x mp x n \X + 1 As x+1 fk x+1 )b x+1 ] (31) 

3 

+ J2l(^s X mp X m \\ + lAs^np^a^ 1 + (XAs x mp x m \X + \As x+1 np x+l )p x+1 ]. 

n 

The characteristic of the two above equations is that only one specific mode in Ath segment (labeling as I or 
to) is left. And the corresponding coefficients is expressed as a function of all the coefficients in (A + l)th 
segment. In addition, if we apply the above procedure to the (A — l)th interface, we can obtain two analogous 
equations. 

a x + e' 2W "e^ k 'y X "b x = ]T[(AB S A Zfc A |A - LBa*" 1 ^- 1 )^" 1 + (XBs x lk x \X - IBs^jk^ 1 )^- 1 } (32) 

3 

+ Y^[(XBs x lk x \X-lBs x - 1 np x - 1 )a x - 1 + (XBs x lk x \X — lBs x ~ 1 np x ~ 1 ) fi^ 1 }, 

n 

and 

Xle^^ai + = "£[(\Bs x m P i |A - lBs^jkf- 1 )^- 1 + {XBs x mp x m \X - \l!, x 'jk x : b x 1 (33) 

3 

+ ^[<AB,s A mp A |A - lBs^np^atT 1 + (XBs x mp x n \X - IBs^np^ 1 )^ 1 ]. 

n 

Combining Eqs. (|30H33p and by a straightforward derivation, it is possible for us to establish an iteration 
relation between the expansion coefficients of the adjacent segments. It is written in a matrix form as 

[C A ] = [M^_ 1 ][C X - 1 ] + [M A +1 ][C A+1 ], 2<\<N-1, (34) 

where [C A ] is defined as a column matrix consisting of all expansion coefficients in Ath segment, namely, 
[C A ] = [a x ■ ■ ■ a x ■ ■ ■ b x ■ ■ ■ a x n ■ ■ ■ {3 x n ■ ■ ■ ] T . The transfer matrices [M^_J and [M A +1 ] are associated with 
the inner products appearing in Eqs. pOH33"]) . We ignore the details about their definition. But it is not a 
difficult task to obtain them from the above equations. The feature of the relation given in Eq. ([5^|) is that 
the coefficient [C A ] is expressed explicitly in terms of the two adjacent ones. If we want to work out the next 
relation, i.e. to express [C A+1 ] as a function of [C A ] and [C A+2 ], we need to use the above procedure at the 
(A + l)th segment. For example, to establish a relation between [C A+1 ] and [C A ] similar to the ones given 
in Eqs. (|32l33p . we need to use the connection condition given by Eq. (f2"3"|) . 
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What we must emphasize herein is that we can not obtain a simpler transfer matrix to connect two sets 
of expansion coefficients between the adjacent segments, namely, in the form of [C x ] — [M^ +1 ][C X+1 ]. This 
is due to that the numbers of modes in the neighbor segments are different. Instead we can only obtain 
an iteration relation involving three adjacent segments, as done above. Just due to such an encumbrance, 
we can not express [C 1 ] in terms of [C 2 ] explicitly. However, by means of the wavefunction connection 
conditions at the first junction interface, we can obtain the following relation 

[Ml][C l ] = [M$)[C% (35) 

Noting that the dimensions of the transfer matrices [Ml] and [M 1 ] are, respectively, J 1 x 2J , and J 1 x 2J 2 . 
Both are not square matrices. By the same token, at the last junction interface we obtain 

[MZ}[C N } = [M»_ 1 }[C»-i]. (36) 

From the above iteration relations, Eqs. (|34j|36p . and by a straightforward derivation we can eliminate all the 
intermediate coefficients, as well as [a 1 ] and [/3 N ]. As a result, we obtain the scattering matrix [S] defined 
inEq.©. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



After establishing the transfer matrix theory to describe the electronic transport through the GMJ, we are 
now in a position to perform the numerical calculation about the conductance spectrum in some typical GMJ 
structures. At first, we consider a simplest case, a graphene single junction formed by interconnecting the 
left and right leads directly. We will focus mainly on the electronic transport property in the vicinity of the 
Dirac point. Therefore, we choose the two leads to be metallic ones which do not present any finite bandgap 
between the conduction and valence bands. In such a case the linear subband across the Dirac point provides 
a basic channel for the electron transmission through the junction. The calculated conductance spectra(<? 
vs E) for two single junction structures are shown in Fig. 2 (a). We see that the spectrum calculated by 
our transfer matrix approach agrees well with the one obtained by the tight binding calculation in the low- 
energy region. In particular, for the junction consisting of relatively wider ribbons the agreement between 
the results of the two kinds of methods is very satisfactory. Basically, this indicates that our generalization 
of the transfer matrix method is successful to describe the electron transport through the graphene junction. 

All these conductance spectra shown in Fig. 2(a) exhibit the staircase-like structures, which can be readily 
explained by the match of the subband structures of the two metallic leads. An interesting point to note in 
Fig. 2(a) is that the conductance spectrum shows a notable suppression in the vicinity of the Dirac energy. 
In particular, a zero conductance occurs at the Dirac energy. At first sight, such a conductance suppression 
is inconsistent with the gapless subband structure of the two metallic leads. The linear subbands of the two 
metallic leads always match each other to provide an electron transmission mode. As a result, one can expect 
that a nonzero conductance appears at the Dirac point. But such a reasoning is dramatically in contradiction 
with the calculated zero conductance. In fact, the nature of the conductance suppression is the antiresonance 
effect, which we can explain as follows. In the vicinity of the Dirac energy, only the linear subbands are 
relevant to the electron transmission. Hence the two metallic GNRs can be viewed as single mode quantum 
wires coupling to each other directly. However, the zigzag-edged shoulder induces a localized edge state 
with the eigen-energy equal to the Dirac energy. Such a localized state couples to the linear subband of the 
wider GNR. Consequently, when the electronic transport is limited in the vicinity of the Dirac point, the 
graphene junction is equivalent to the T-shaped quantum dot structure as shown in Fig. 2(b). The linear 
conductance of such a model has been extensively studied ) 48 ! 49 ' 50 which can be expressed in terms of the 
relevant parameters 

G(E) - ^ T{E) ~ T(i + iF( E -er + [^r (37) 

where T = 2irp2T 2 and £ = it 2 p\piv 2 with P\(2) being the electron density of the states in two leads. This 
expression presents a zero conductance at the quantum dot level e indeed, which is called the antiresonance 
effect. The antiresonance is in fact a destructive quantum interference. The lateral quantum dot introduces 
new Feynman paths with a phase shift tt. As a result, the destructive quantum interference occurs among 
electron Feynman paths. In the metallic graphene junction, the edge state attached to the zigzag-edged 
shoulder of the junction plays the role of the laterally coupled quantum dot, which results in the antiresonance 
at the Dirac point. In Fig. 2(c) we calculate several conductance spectra when the width of the narrower lead 
approaches the wider one whose width is fixed. We find that the suppressed conductance valley becomes 
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narrower as the widths of the two leads approach each other. This effect can be explained in terms of 
the above quantum dot model. One can readily infer that the parameters v(the coupling between the two 
leads) becomes larger, but 7"(the coupling between the propagating mode in the wider lead and the edge 
state localized at the junction shoulder) becomes smaller, as the two leads approach each other in their 
widths. Consequently, from Eq. (|37[) one can deduce that the conductance valley becomes narrower. On the 
contrary, when we fix the width of the narrower lead, but to increase the width of the wider lead, we can 
see from Fig. 2(d) that conductance valley does not vary notably. This is because that in such a situation 
the parameter r, which mainly controls the valley width, does not vary notably as the wider lead is further 
widened. 

In Fig. 3(a) we show the calculated conductance spectrum for a step-shaped GMJ. We take the case of 
N = 3, namely, a two-junction structure, slightly more complicated than a graphene single junction. First, we 
can see that the results calculated by the tight binding and transfer matrix methods, respectively, agrees with 
each other very well. But the tight binding calculation is much more time-consuming than the corresponding 
calculation via the transfer matrix method. Secondly, we can sec that the conductance valley around the 
Dirac point found in a graphene single junction remains. And it is further broadened. In Fig. 3(b), the width 
of the right lead(the wider one) is fixed, the calculated spectrum shows that the conductance valley gets 
narrower as the two shoulders becomes shorter. In contrast, if we fix the width of the left lead, meanwhile, 
to increase the width of the two shoulders, the conductance valley is almost invariable. Such a case is 
shown in Fig. 3(c). These results shown in Fig. 3(b) and (c) is consistent with the size dependence of the 
graphene single junction, as shown in Fig. 2(c) and (d). This indicates that the conductance characteristic 
in a single junction is preserved and further strengthened in the two-junction structures. The results of the 
conductance spectra shown in Fig. 3(d) indicate that the profile of the steep conductance valley does not 
vary sensitively with the increase of the longitudinal size of the intermediate ribbon segment. This reflects 
that the interaction between the two edge states at the two adjacent shoulders is very weak. 

The conductance spectra for more complicated step-shaped and T-shaped GMJs are shown in Fig. 4(a) 
and (b), respectively. A common feature of these spectra is that the conductance valley gets steeper as the 
number of the junctions increases. As shown in Fig. 4(a) and (b), when N > 3 the profile of the conductance 
valley becomes stable, unsensitive to the further increase of the number of the junctions in the structures. 
This indicates that the conductance valley in a single junction develops rapidly into a well-defined insulating 
band in the GMJs. Following the antiresonance mechanism mentioned above, it is not very difficult to 
understand the appearance of the well-defined insulating band in the GMJ. As claimed above, the graphene 
single junction in the low-energy electron transport regime is equivalent to a lateral quantum dot structure. 
In an analogous way a mapping between the GMJ and a laterally coupled quantum dot chain, which is 
illustrated in Fig. 4(c), is reasonable. Of course, the edge states localized at the individual shoulders of 
the junction interfaces play the roles of the dangling quantum dots. And the quantum dots in the main 
chain mimic the electron states in the ribbon segments. Each quantum dot pair, made of the quantum dot 
in the main chain and its lateral attachment, can be viewed as an antiresonance unit which contributes a 
transmission probability amplitude, the squared modulus of which has the same form as the transmission 
probability given by Eq. (j3"T|) . In the GMJ the electron wave will be reflected many times by the multiple 
interfaces before it finally enters into the right leads. The shortest transmission path is the one that the 
electron penetrates all the junctions directly without any reflection. Even in such a situation, the electron will 
undergo N junctions, in other words, N antiresonance units. Accordingly, the corresponding transmission 
probability is proportional to the iVth power of T(E), i.e. T(E) N . (T(E) is the transmission probability 
of an individual antiresonance unit given by Eg. (1371) V Thus the conductance suppression around the Dirac 
point included in T(uj) of each quantum dot is efficiently strengthened. As a result, the antiresonant valley in 
a single junction evolves into a well-defined insulating band in a GMJ due to multiple quantum interference. 
A comprehensive explanation about the occurrence of the insulating band in a lateral quantum dot chain 
was presented in the language of Feynman path in our previous work ) 50 i 51 which can be reasonably used to 
account for the well-defined insulating band in the present GMJ structures. 

In Fig. 5 we show the conductance spectrum which is adjusted by applying a gate voltage under some 
junctions interfaces. In our theoretical treatment, the gate voltage is mimicked by introducing a finite 
potential constant in the ribbon segment where the gate voltage is applied. We can see that the insulating 
band is effectively broadened when certain of the junction interfaces are covered by the gate voltage. But 
the cost of such a broadening is that the insulating band is not so well-defined than the one in the absence 
of a gate voltage. Such a result can be readily understood according to the multiple quantum interference 
picture explained above. The transmission probability of each junction possesses a valley center around the 
level of the localized edge state, in other words, the quantum dot level in the corresponding quantum dot 
chain. In the absence of a gate voltage such a level is just at the Dirac point. On the other hand, in the 
presence of a gate voltage covering partially some junctions of the total structure, the levels of the edge 
states of these junctions, namely, the centers of these antiresonance units, separates from those of other 
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junctions free from the gate voltage. Accordingly, the total transmission probability presents a broadening 
of the insulating band, due to that the gate voltage destroys the superposition of the antiresonance centers 
of all antiresonance units. By the same argument, one can readily understand that the insulating band only 
exhibits a shift without an observable broadening when all junctions are simultaneously tuned by a gate 
voltage. In light of the sensitive adjustment of a gate voltage on the insulating band, as shown in Fig. 5, we 
suggest that the GMJ can be considered as a device prototype of a nanoswitch. 

IV. SUMMARY 

By generalizing the transfer matrix approach which was established to describe the electron transport 
through quantum waveguide made of conventional semiconductor materials, we provide an alternative theo- 
retical method to calculate the conductance spectrum of a GMJ, in place of the tight-binding treatment. By 
comparing the calculated results by the two methods, we find that our generalized transfer matrix method 
works well in describing the electronic transport properties of the GMJ structure if its size is not very small. 
In particular, this method is especially suitable to deal with the GMJs with relatively large sizes and many 
junctions. In such a case, it is much more timesaving to calculate the conductance spectrum by this method 
than by the tight-binding calculation. Interestingly, with the transfer matrix method, we find that the GMJ 
exhibits a well-defined insulating band around the Dirac point, which is due to the strengthen of the intrinsic 
antiresonance effect of a graphene single junction by the multiple quantum interference at many junction 
interfaces. By virtue of the feature that the insulating band can be sensitively tuned by a gate voltage, we 
propose to design a nanoswitch based on this GMJ structure. Our generalized transfer matrix method is also 
suitable to treat the graphene waveguides with various junction shapes, such as the tortuous nanoribbons 
with distinct bent angles. 

Finally, we need to make some remarks on the issue to implement the function of nanoswitch, by means 
of the tunability of the insulating band of the GMJ by a gate voltage. In fact, such a transport phenomenon 
also appear in the T-shaped waveguides and the laterally coupled quantum dot chainS ) 51 i 52 made of the 
conventional semiconductor materials. The similar proposal to make a nanoswitch was previously mentioned 
in these relevant work. However, the occurrence of a well-defined insulating band in these quantum structures 
requires that the fluctuation of the quantum dot energy levels (the size of the stubs in a T-shaped waveguide) 
to be sufficiently small. It is indeed a challenging task to fabricate many identical quantum dots and 
connecting them periodically in a circuit to realize the function of nanoswitch. However, in the present GMJ 
structures, the problem of energy level fluctuation is automatically avoided since the antiresonant levels are 
provided by the localized edge state at the zigzag edge, rather than the quantum confinement in the stub 
region. The energy levels of these localized states are automatically aligned with each other at the Dirac 
point, regardless of the details of the shape of the GMJ. This implies the prefect realization of the uniform 
antiresonance levels. Generally, the insulating band can be easily formed in the step-shaped as well as the 
T-shaped GMJs. It does not rigidly require the identity of the shapes of every junctions. In addition, it 
is experimentally demonstrated that ballistic transport can be retained in graphene nanostructures over 
sub-micron scale. This means that the quantum interference can play a dominate role even the size of the 
GMJ becomes very large. Thus the insulating band structure is expected to remain in a GMJ even with 
a size up to sub-micron scale. To sum up these feature, we can conclude that the GMJs are the optimal 
candidate as the prototype of the nanoswitch. 

This work was financially supported by the National Nature Science Foundation of China under 
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FIG. 1: (Color online) (a) Schematic illustration of a graphene multiple junction structure, (b) The honeycomb lattice 
of an intermediate part of the graphene multiple junction, around the Ath junction. The carbon atoms belonging to 
the two distinct sublattices, A and B, are distinguishingly labeled as A:» and B:o. The atoms at the edges linked to 
the interior ones by the dashed lines stand for the truncated atoms to shape the GMJ. 

FIG. 2: (Color online) (a) The conductance spectrum of two graphene single junctions with distinct sizes. Comparison 
is made between the conductance spectra calculated by the tight binding(TB) and the transfer matrix(TM) methods, 
(b) Schematic of a laterally coupled quantum dot structure, where v and r represent, respectively, the coupling 
coefficients between the two leads and between the quantum dot and the right lead. The calculated conductance 
spectrum from this model is also plotted, (c) A comparison of the conductance spectra of the graphene single 
junctions by changing the width of the left lead, whereas the width of the right lead is fixed, (d) Similar to (c), but 
the width of the right lead is changed with the left lead fixed. 

FIG. 3: (Color online) Conductance spectrum of step-shaped two-junction structures(GMJs of N=3). (a) A compar- 
ison of the calculated results by the tight binding and the transfer matrix methods, (b) and (c) A comparison of the 
conductance spectrum by changing the width of the junction shoulders. The two different cases with the width of 
the left and right leads fixed, respectively, are shown in (b) and (c). (d) A comparison of the conductance spectrum 
by changing the longitudinal size L2 of the intermediate segment. 

FIG. 4: (Color online) Conductance spectrum of typical GMJ structures. The length of the intermediate ribbon 
segment is L\ — 48/\/3. (a) The result of the step-shaped structures, (b) The result of the T-shaped structures, (c) 
A illustration of the laterally coupled quantum dot chain, which is used to account for the insulating band in the 
conductance spectrum of the GMJ. 



FIG. 5: (Color online) Conductance spectrum of GMJ structures in the presence of the gate voltage, (a) The result 
of the step-shaped structures, (b) The result of the T-shaped structures. V\ denotes the potential constant produced 
by the gate voltage in the A segment. 
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